Triplet supercurrent due to spin-active zones in a Josephson junction 
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Motivated by a recent experiment evidencing triplet superconductivity in a ferromagnetic Josephson junction 
with a Cu2MnAl-Heusler barrier, we construct a theoretical model accounting for this observation. The key in- 
gredients in our model which generate the triplet supercurrent are spin-active zones, characterised by an effective 
canted interface magnetic moment. Using a numerical solution of the quasiclassical equations of superconduc- 
tivity with spin-active boundary conditions, we find qualitatively very good agreement with the experimentally 
observed supercurrent. Further experimental implications of the spin-active zones are discussed. 
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Since the pioneering studies 50 years ago 1 , the interplay 
between ferromagnetic (F) and superconducting (S) order has 
been much investigated, particularly so in recent yeard^E. This 
can largely be ascribed to important experimental develop- 
ments which have allowed for microscopic studies of both 
artificially engineered 4 and intrinsic coexistence^ of F and S 
order, in addition to theoretical advances. One of the most ex- 
citing prospects in hybrid F|S structures is the possibility of 
tailoring the desired properties of the system on a nanometer 
scale. To accomplish this, it is necessary to take seriously the 
influence of the interface properties. Depending on whether 
the interfaces have spin-dependent properties or not, exotic 
new features may come into play in F|S structures, including 
long-range Josephson effects^' and unconventional types of 
superconducting pairing 9 *^. 

In a very recent experiment by Sprungmann et al!^\ the 
importance of such interface properties was underscored. 
A long-range triplet supercurrent was observed in a S|F|S 
Josephson junction with a Cu2MnAl-Heusler barrier acting 
as the F region. Surprisingly, Sprungmann et al. observed a 
conventional, exponentially decaying supercurrent in the as- 
prepared state of the Cu2MnAl layer, whereas a long-range 
supercurrent virtually independent of the junction thickness 
(up to a critical value) was observed upon annealing the junc- 
tion. Beyond a critical value of the junction-thickness, an 
unusual abrupt decay of the supercurrent was observed. In 
the annealed case, the Heusler layer acquires a ferromagnetic 
order in its core, but retains spin glass order near the inter- 
faces in the thickness range just above the onset of ferromag- 
netism. Therefore, it was suggested in RefP 

that the cou- 
pling between ferromagnetic and spin glass order would lead 
to a canted magnetization texture near the interfaces, which 
could be the necessary mechanism responsible for the ob- 
served triplet supercurrent. Illuminating this matter would be 
important to understand further the role of spin-active zones 
in Josephson junctions and their possible manipulation, which 
in turn could lead to tunable long-range supercurrents. 

In this Rapid Communication, we construct a theoretical 
model to explain the experimental finding in RefPby includ- 
ing the role of canted magnetic interface moments in the qua- 
siclassical theory of superconductivity. Employing a numer- 
ical solution in the diffusive regime of transport, we obtain 
very good qualitative agreement with observed supercurrent 
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FIG. 1: (Color online) The model employed for a Josephson junction 
with a ferromagnetic Heusler Cu2MnAl barrier. The junction width 
is L, and we take into account a canted magnetization texture near the 
interfaces with misalignment angles o.l,r relative the bulk magneti- 
zation. These spin-active zones generate a long-range supercurrent. 



in RefP. Moreover, we highlight additional implications of 
spin-active zones in a Josephson junction. In particular, we 
find that if these zones couple together in a parallell align- 
ment, the total supercurrent may actually vanish for certain 
misalignment angles of the canted magnetic moments. Our re- 
sults demonstrate the significance of spin-active zones in fer- 
romagnetic Josephson junctions, and may provide a guideline 
for interpretation of experimental data and future investiga- 
tions of such heterostructures. 

The system under consideration is shown in Fig. [T] A ferro- 
magnetic region of width L is sandwiched between two stan- 
dard i-wave superconductors (e.g. Nb). Near the interfaces, 
we identify spin-active zones which may arise e.g. due to the 
presence of magnetic disorder or canted magnetic moments 
which are misaligned compared to the bulk magnetization di- 
rection. To study this system, we employ the quasiclassical 
theory of superconductivity which provides equations of mo- 
tion for the Green's functions of the system. By supplement- 
ing these equations with proper boundary conditions for spin- 
active interfaces with possibly misaligned magnetic moments, 
we are able to compute the supercurrent flowing through the 
system. We give a brief account of the theoretical framework 
here, and refer the reader to e.g. Ref. 2 for a comprehensive re- 
view. In the F region, we employ a Ricatti-parametrizatiorP 
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of the Green's function as follows 
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where 9\[ = (I + yy) 1 and 3\£ = (1 + yy) 1 . It also satisfies 

(g) 2 = 1, and is characterized by the two unknown 2x2 ma- 
trices y and y. The above Green's functions satisfy the Usadel 
equatiorJEl 



£>gd x g + /[ep 3 + Mag (x 3 , x 3 ) , |] = 0. 



(2) 



Here, T> is the diffusion coefficient in the F region, h is the 
magnitude of the exchange field, while e is the quasiparticle 
energy. To obtain a complete solution of the Green's function, 
it is necessary to supplement the Usadel equation with bound- 
ary conditions^. Assuming that the F region extends from 
x = to x = L, the boundary conditions read 

2i T Lgd x g = [§',§} +it\fi(a L ),g] ztx = Q, 
2f T Lgd x g = [g, g r ] - if s [M(a R ) ,g]atx = L. (3) 



Here, we have introduced 

M(a) = (cosa)diag(a z ,aJ + (sina)diag(a v ,a v *) 



(4) 



with (j) being the angle between the barrier magnetic moment 
and the bulk magnetization (z-axis) as indicated in Fig. [T] 
Above, Oj denotes the j'th Pauli matrix in spin space. We 

have defined "/j r = R l g /Rf and "// = 1/ (R^G 1 /), where Rn 

l r 

is the normal-state resistance of the N region, R B ' is the bar- 

/ r 

rier resistance at the left/right interface, while G s ' accounts 
for the spin-dependent interfacial phase-shifts at the left/right 
interface. In the S region, we make use of the bulk solution 
for the Green's function 
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where c = cosh 9, s = sinhG, 8 = atanh(|Ao|/e), and the super- 
conducting phase difference is %■ This approximation is valid 
for low transparency interfaces in which case the proximity 
effect is weak. The above equations constitute a closed set 
which may be solved numerically, and we add a small imag- 
inary part to the quasiparticle energy for improved numeri- 
cal stability, i.e. e — > e + i8 with 8/Ao = 0.01. To model 
the experimental situation in Refill, vve consider the weak 
proximity effect regime with a low-barrier transparency and 
an exchange field h. Specifically, we choose Yy' = 40 and 
h/Ao = 80 which is within the regime of validity for qua- 
siclassical theory when assuming Ao ~ 1 meV. Also, in this 
regime one may safely neglect magnetoresistance terms Jmr 
in the boundary conditions 16 . We have chosen a numerical 
approach here to avoid overburdening the paper with cumber- 
some analytical expressions. 



It was proposed by Sprungmann et al. ^ that the micro- 
scopic mechanism generating the triplet supercurrent was the 
specific magnetization profile in the Cu2MnAl layer, featuring 
a canted magnetization texture near the interface regions. An 
important observation in RefP^ was that the triplet supercur- 
rent was only observed in a limited thickness range, and died 
off rapidly above a critical thickness L c thus leaving behind 
only the conventional short-ranged Josephson current. This 
indicates that the misalignment angle of the interface mag- 
netization compared to the bulk depends on the thickness L 
of the layer. Such a conjecture agrees with the fact that the 
Heusler layers gradually make a transition near the interfaces 
from pure spin-glass to coexisting spin-glass order with a fer- 
romagnetic moment of increasing magnitude as the thickness 
increases. To phenomenologically model such behavior, we 
write the misalignment angles shown in Fig. [T]as 



a = ao 



(6) 



where ao is the misalignment angle at the onset of ferromag- 
netic order and K, 1 accounts for the slope with which the 
canted moments relax into the same orientation as the bulk 
magnetization. 

Let us further qualify 17 the model employed here for the 
canted interface magnetization texture. A thin film Heusler 
alloy features a spin glass structure if the film is sputtered at 
room temperature. Annealing will induce a transposition of 
nearest-neighbour atoms and tries to drive the compound into 
a ferromagnetically coupled structure. The maximum order- 
ing effect is achieved within the core of the layer, whereas 
at the interfaces this ordering process is disturbed due to e.g. 
interdiffusion. Here, an antiferromagnetic Mn-Mn-coupling 
will persist locally and competes with the ferromagnetic or- 
dering. This causes the mentioned canting of the magnetiza- 
tion axis defined by the one of the core. The issue of what 
the magnetization profile looks like at the interfaces of thin 
Heusler layers was investigated experimentally in Ref!^. We 
here also briefly mention a theoretical work on the influence 
of disorder in Co-based Heusler alloys by Picozzi et al\^\ al- 
though the comparison should be considered carefully since 
in Co-based Heuslers the Co contributes effectively to ferro- 
magnetism whereas Cu does not carry any magnetic moment 
in the relevant structure. 

We are now in a position to calculate the normalized critical 
current density of the junction, which in the weak-proximity 
effect regime is obtained via 
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where we have defined 

A M ± A = [/±(-e)]*3,/ T (e) -/±(e)3,[/ T (-e)]*, 
AoM CT /£ = [/b(-e)]*dx/«(e) - Me)d x [f c (-e)}* . (8) 

Here, {/f,/i,/-j-,/|} denote the S z = triplet, singlet, and 
equal spin-pairing anomalous Green's function induced in the 
F region, and we have defined f± = f,±f s . To make direct 
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FIG. 2: (Color online) Plot of the critical current density vs. the 
junction width L. Triangles (black) represent results for the case a = 
0, i.e. no canted interface moments. Squares (red) represent results 
for misalignment angle a given by Eq. [6] We have set ai = oc« 
with oio = Jl/4 and ^/ = 15. Assuming a superconducting coherence 

length of ^ ~ 14 nm in dirty Nb'21 we use Lrji, = 0.76 to obtain 
the sharp drop-off near L = 10.5 nm as in RefU^I Inset: Zoom-out 
version of the current density. 



contact with the experiment of RefP^l we contrast a situation 
with a = 0, i.e. no canted interface moments, against a situ- 
ation where the misalignment angle a evolves with the thick- 
ness L according to Eq. Q. The result is shown in Fig. [2] The 
result is seen to be qualitatively in very good agreement with 
the finding or Ref. 12 : the supercurrent shows little decay upon 
increasing the thickness L in the regime where the spin-active 
zones generate a misalignment angle, and then collapses onto 
the conventional singlet result above a certain thickness L c . 
This suggests that the spin-dependent interface properties may 
play a pivotal role in the generation of the triplet supercurrent. 
In our picture, the abrupt change in the spin-triplet Josephson 
current is due to an abrupt change in the angle a in Eq. [6] 

Another interesting finding in Ref. 12 is the non-monotonic 
temperature dependence of the critical current observed in the 
region of widths L where the long-range current is dominant. 
With T c ~ 8 K being the critical temperature 17 , we plot the 
current vs. temperature in Fig. [3] for three choices of the 
width: (i) right below the region of a dominant triplet current 
(L = 6.1 nm), (ii) in the middle of the region of a dominant 
triplet current (L = 8.1 nm and 8.6 nm), and (iii) right after 
the vanishing long-range current (L — 10.7 nm). We again 
obtain a good match to the results of Ref. 12 . It should be men- 
tioned that we do find signs of 0-TC oscillations for some partic- 
ular choices of intermediate widths L, but the non-monotonic 
behavior in the regime of a dominant long-range current is 
nevertheless confirmed. Close examination of the current vs. 
temperature curves in Ref. 12 reveals a hint of small-scale os- 
cillation superimposed on the overall non-monotonic behavior 
for some particular widths L. We were not able to reproduce 



FIG. 3: (Color online) Plot of the critical current density vs. temper- 
ature. To make contact with the experiment in Ref. 12 where T c ~ 8 
K, we normalize the current to its value at T /T c = 0.25. The rest of 
the parameters are as in Fig. [2] 
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FIG. 4: (Color online) Plot of the critical current density vs. 
the misalignment angle a in the spin-active zones. In (a), 
we set O.L = 0C/e while in (b) (Xi = —Ur. From top to bot- 
tom at a = 7t/2, the curves correspond to thicknesses df /i; = 
{0.60, 0.62, 0.64, 0.66, 0.68, 0.70}. 



these fine-scale oscillations within our model, and speculate 
that these might pertain to a more complicated magnetization 
profile in the bulk of the Cu2MnAl layer since similar behav- 
ior has been predicted 20 in the conical ferromagnet Ho. We 
intend to address this issue in a forthcoming work. It should 
be noted that the spin-dependent phase-shifts occuring at the 
interfaces can in general influence the 0-71 transition pattern 
both as a function of the width and temperature of magnetic 
as well as non-magnetic (normal interlayer) Josephson junc- 
tions. 

To further highlight the influence on the current density by 
the presence of a canted magnetization texture near the inter- 
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faces, we plot in Fig. |4]the current density vs. the misalign- 
ment angle in the case of an (a) parallel coupling (Xl = Or 
between the spin-active zones and an (b) antiparallel coupling 
Of, = — (Xr. The main difference between these two cases is 
that the total supercurrent may actually vanish at some mis- 
alignment angles when the coupling is parallel. In this sce- 
nario, the contribution from the long-range triplet current ef- 
fectively cancels out the singlet supercurrent. Common for 
both scenarios is that the maximally attainable critical current 
occurs near a = Jl/2. However, it should be noted that the 
maximum does not occur exactly at a = n/2, but rather at a 
slight off-set a c < n/2. When we reverse the bulk magnetiza- 
tion direction, i.e. h — > (—h), we note that the maximum of the 
critical current occurs at an angle a c > n/2 which is equidis- 
tant from n/2 compared to the former case. This behavior can 
be understood from a symmetry perspective. Namely, the sys- 
tem is not invariant under a spatial inversion whenever h ^ 
due to the definite direction of the magnetization. Therefore, 
the angle a c providing the maximum critical current will ei- 
ther be smaller or larger than n/2, depending on whether h 
points along z or — z. 

So far, we have constructed a theoretical model for the spin- 
active zones in the Heusler layer which is able to account well 
for the experimental finding in Refill There are nevertheless 
extensions of this model which could be appropriate to pursue 
in order to further clarify the underlying physics, and which 
we comment on. In our treatment, we have included the spin- 
active zones as effective canted magnetization textures near 
the interface and included the resulting conductance-terms in 
the quasiclassical boundary conditions. This is certainly only 
an effective model for the interplay between spin-glass and 
ferromagnetic order in the Heusler layers in the interface re- 



gion. In a more microscopic lattice-model treatment of the 
system, one could address how the specific details of such an 
interplay would influence the critical current behavior. An- 
other issue of interest, is the precise role of the spatially de- 
pendent magnetization texture in the ferromagnetic layer. Re- 
cent theoretical works 21 have highlighted the necessary re- 
quirements of the magnetization profile which renders possi- 
ble a long-range triplet current. It would thus be interesting to 
experimentally extract this profile and its dependence on both 
spatial coordinate and its evolution in terms of magnetization 
direction, so that one could construct a more accurate model 
for this system. Qualitatively, we expect that our model cap- 
tures the essential features which could account for the gener- 
ation of a long-range triplet supercurrent in the experimental 
system of RefM^l 

In summary, we have constructed a theoretical model to ac- 
count for recent experimental findings in RefP which pre- 
sented evidence for triplet superconductivity in a ferromag- 
netic Josephson junction with a Cu2MnAl-Heusler barrier. 
The crucial ingredients in our model which generate the triplet 
supercurrent are spin-active zones near the interfaces due to a 
coupling between ferromagnetic and spin glass order, result- 
ing in an effective canted magnetization texture near the inter- 
faces. Using a numerical solution of the quasiclassical equa- 
tions of superconductivity with spin-active boundary condi- 
tions, we find qualitatively very good agreement with the ex- 
perimentally observed supercurrent. Further experimental im- 
plications of the spin-active zones have also been discussed. 
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